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Abstract 

In this paper, the physical realizability condition 
for a specific class of nonlinear quantum systems is re- 
lated to the lossless property of nonlinear dissipative 
systems having a specific storage function. 



1. Introduction 

There are two types of quantum stochastic differ- 
ential equations, namely, linear and nonlinear. Linear 
quantum systems can be described by linear quantum 
stochastic differential equations and arise mostly in the 
area of quantum optics; e.g., see Q, |9) and (SJ. An 
important class of linear quantum stochastic models de- 
scribe the Heisenberg evolution of the annihilation and 
creation operators of several independent open quantum 
harmonic oscillators that are coupled to external coher- 
ent bosonic fields, such as coherent laser beams; e.g., 
see (8], [7 1 and [9 |. A special class of these linear quan- 
tum systems is driven by quantum Wiener processes as 
in where physical realizability conditions are devel- 
oped to determine when the linear quantum system un- 
der consideration can be regarded as a representation 
of a linear quantum harmonic oscillator. This imposes 
some restrictions on the matrices describing the linear 
quantum model. This notion of physical realizability 
has been further investigated in ifTTIl where the authors 
prove the equivalence between the algebraic conditions 
for physical realizability obtained in [ 2 1 and a frequency 
domain condition that an associated linear system is 
(7,7)-unitary. In fact, ifTTl extends the frequency do- 
main physical realizability results of fl2l . iPTSIl and JT] 
which apply to linear quantum systems described purely 
in terms of annihilation operators. More explicitly, in 
fl~2), IIT31 and (TJ, the physical realizability conditions 
developed for annihilation-operator linear quantum sys- 



tems have been related to the lossless property of lin- 
ear systems and [11] generalizes this result to relate the 
corresponding physical realizability conditions of anni- 
hilation and creation operators linear quantum systems 
to the (7,7)-unitary property. 

In this paper, we restrict attention to a class of non- 
linear quantum systems which is a generalization of the 
annihilation operator only linear quantum systems con- 
sidered in [12 1, [13 1 and [1[. Then, we relate the re- 
sulting physical realizability conditions to the lossless 
property of nonlinear dissipative systems having a spe- 
cific type of storage function. 

2. The Class of Nonlinear Quantum Sys- 
tems 

We consider an open quantum system G with phys- 
ical variable space s^q consisting of operators defined 
on an underlying Hilbert space Mq. The self-energy 
of this system is described by a Hamiltonian 3%? £ s^g- 
The system is driven by m field channels given by the 
quantum stochastic process W 



W 



w„, 



These describe the annihilation of photons in the field 
channels and are operators on a Hilbert space F, with 
associated variable space '. In that case, F is the 
Hilbert space defining an indefinite number of quanta 
(called a Fock space |5|), and & is the space of opera- 
tors over this space. 

We assume the process W is canonical, meaning 
that we have the following second order Ito products: 
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dW k (t)dW,(t)* 
dW k (t)*dWi(t) 
dW k {t)dWi(t) 
dW k (t)*dWi(t)* 



0; 
0; 




where Wk(t)* is the operator adjoint of W J t(f) which is 
defined on the same Fock space. 

The system is coupled to the field through a scat- 
tering matrix S = I and a coupling vector of operators L 
given by 



L = 



where Lj £ s^q. 

The notation G = (S,L,H) is used to indicate an 
open quantum system specified by the parameters S,L 
and H where H represents the Hamiltonian of the sys- 
tem. For S = I, the Schrodinger equation is given by 

dU(t) = {dW Jt L-LUW-~L t Ldt- Mdt^U(t) 

(1) 

with the initial condition U (0) = /. Equation (Q]i deter- 
mines the unitary motion of the system in accordance of 
the fundamental laws of quantum mechanics. Note that 
the notation ^ refers to the Hilbert space adjoint. 

Given a system annihilation operator a/ S s^q, its 
Heisenberg evolution is defined by a/(f) = j t {ai) = 
U (t)*aiU (t) and satisfies 

da,(t) = (& L (t)(ai(t))-i[ai(t),H(t)]) ( 2 ) 
+ [L(t)\a,(t)] dW(t) 

where the notation [A,B] = AB — BA denotes the com- 
mutator of two operators A and B. 

In equation (0, all operators evolve unitarily and 
the notation Jz?i(a/) refers to 

& L (ai) = ^tf[ai,L] + ^[tf,ai\L. (3) 

In this paper, we restrict attention to annihilation only 
coupling operators; i.e, the coupling operator L{t) de- 
pends only on the annihilation operator a/(f) and not 
the creation operator a/ (f)*. Therefore, L(t) satisfies the 
following commutation property: [a/,L] = and then 

J? L (a l ) = ^[lS,a l ]L. (4) 

The generator of the system G is given by 

<g G {a l ) = -i[a u H]+& L {a l ). (5) 

For the case of having n annihilation operators 
a\,... ,a n , we define 



& u {a) = -[L\a l ]L 



and 



& G .(a) = -i[a l ,H] + £' L .(a). 
Therefore, we can write 

dai{t) = {& Li {t)(ai{t))-i[ai(t),H{t)\ 
+ [L{t)\a t {t)]dW{t) 



(6) 
(7) 



and 



da(t) = (& L {t){a(t))-i[a(t),H{t)\ 
+ [L(t)\a(t)] dW(t) 



where 



and 



-%(a) = 



[a,H] = 




(8) 



(9) 



Let 



A(a(f),a(f) t ) 



/ A 1 (a(f),fl(f) t ) 

V A„( fl (f), fl (f) t ) 
2iM 



= S£ L (a)-i[a,H] 




and 



B(a(t),a(rf) 



/ Bi(a(f),a(f) t ) 

V B„(fl(f),a(f) t ) 
/ [L T ,ai] 

~L\a] . 



(10) 



(11) 



(12) 



Hence, 

da(t) 
da(t)* 



A(a(t),a(t) f )dt + B(a(t),a(t) f )dW(t) 
A(a(t),a(t) f )*dt+B(a(t),a(rf)*dW(ty 



Note that for the case of matrices, the notations * and 
t refer respectively to the complex conjugate and the 
complex conjugate transpose of the matrix in question. 

(13) 

On the other hand, the components of the output fields 
are defined by y(t) = j t (W(t)) = U (t)*W (t)U '(f) and 
satisfy the nonlinear quantum stochastic differential 
equations 

dy(t) = C(a(t))dt + D(t)dW(t) 
dy(t)* = C(a(t))*dt+D(t)*dW(t)* 

(14) 

where C{a{t)) = L(t), C{a{t))* = L(t)*, D(t) = I and 
D(t)* = I. Hence, the system G can be described by 
the following nonlinear quantum stochastic differential 
equations 



da(t) = 
dy(t) = 

where a(t) 



A(a(t),a(t) ' )dt+B(a(t),a(rf)dW(t); 
C{a(t),a{t) t )dt+D{a{t),a(t)^)dW({).5) 



a 
a 



A(a,a*) = 



C(o,fl T ) = 



C{a) 
C(a)* 



A{a,cfi) 
A(a,cr) 

B(a,a^) 



B(a,a f ) 
B(a,a f ) 
dW(t) 
dW{t)* 



,D(t) = 



D(t) 
D(0* 



and 



dW{t) 



Definition 2.1 The class of nonlinear quantum sys- 
tem we consider in this paper are nonlinear quantum 
systems that can be represented by the QSDES A15\l 
such that the matrices A{a,a'),B{a,a^) and C(a) sat- 
isfy \B(a,cfi),a T '\ =0, [(^(a),^! 7 "] =0, [A(a,a' = 
— [a,A(a,a^) r ] and 

m kj m hi „ n 

<r,=0/i,=0/=l p=\ 
m k v n 

cm) = E L&vp«p ( 16 > 

k v =0p=l 

for i = 1 , . . . , n and v — 1 , . . . , m with m^., m^., m^, fik v p 
and OLpikjhj integers. 

Moreover, the matrices A{a, a*), B(a,a^) andC(a) 
satisfy the following property: 



1 



A(a,a T )^9 x a,a \a [ 9 'Afa, 

i 2n 

1 



A(a,a^) — — B(a,a^)C(a,fl^) 



(17) 



where n = sup (hi + hi) 

\<i<n 







© 

0* 



and = [a,a^]. For more information, on 0, refer to 
Section III below. 

In the sequel, we assume that dW(t) admits the fol- 
lowing decomposition: 

dW(t) = f} w (t)dt + dw(t) 

where w(t) is the noise part of W(t) and j3 vv (f) is an 
adapted process (see [5 1 and [6|). The noise w(t) is a 
vector of quantum Weiner processes with Ito table 



dw(t )dw(t) Jf = Fadt 



(18) 



(see 0) where F$ is a non-negative definite Hermitian 
matrix. Here, the notation t represents the adjoint trans- 
pose of a vector of operators. Also, we assume the fol- 
lowing commutation relations hold for the noise com- 
ponents: 

[dw{t),dw{tf] =dw{t)dw{rf -{dw{t)*dw(tf) T = T w dt. 

(19) 

Here T w is a Hermitian matrix and the notation T de- 
notes the transpose of a vector or matrix of operators. 
The noise processes can be represented as operators on 
an appropriate Fock space; for more details, see [5]. 

The process j3 w (f) represents variables of other sys- 
tems which may be passed to the system $15[ via an 
interaction. Therefore, it is required that j3 vv (0) be an 
operator on a Hilbert space distinct from that of ao and 
the noise processes. We also assume that j5 w (t) com- 
mutes with a(t) for all t > 0. Moreover, since /3 w (f) is 
an adapted process, we note that j3„, (f) also commutes 
with dw{t) for all t > 0. We have also that [a,dW t ] = 0. 

3. Commutation Relations 

For the nonlinear quantum system ( fT5l ). the initial 
system variables a(0) = ao consist of operators satisfy- 
ing the commutation relations 

[fly(0),aj(0)] = & Jk , M=l, ...,«. 
[aj(0),at(0)] = Q jk , 7,fc=l,...,2n.(20) 

Here, the commutator is defined by [flf,fl£] ~ a j a t ~ 
a\aj = © jk where is a complex matrix with elements 
&jk- With a T = (a i , . . .,a„), the relations (l20b can be 
written as 



tl A 



aa~ — (a* a 

© 

©* 



T\T 







= 0. 



(21) 



3.1. Preservation of the Commutation Rela- 
tions 

The following theorem provides an algebraic char- 
acterization of when the nonlinear quantum system 
< fT3T > preserves the commutation relations (l20b as time 
evolves. 

Theorem 3.1 For the nonlinear quantum system 
( 1751 ), we have that [a p (0),a*(0)] = ® pq implies 
[dp(t),dq(t)] = & P g for all t > with p,q — I . . . ,2n 
if and only if 

[A (a, a 1 "), a 1 "] + [a,A(a,fl t ) t ] 
+B(a,a t )T^B(a,a f ) t =0; 
\S(d,a^),a^] = and 
[a,B(a,a f ) f ] =0 (22) 

•it. T \ T *> 

with T^= n T * ■ 

w 

4. Physical Realizability And The Nonlin- 
ear Quantum Harmonic Oscillator 

In this section, we consider conditions under which 
a nonlinear quantum system of the form (TT3T > corre- 
sponds to a nonlinear open quantum harmonic oscilla- 
tor. Such a system is said to be physically realizable. 
The class of nonlinear open quantum harmonic oscil- 
lators under consideration are defined by the nonlinear 
coupling operator L and a nonlinear Hamiltonian 77. To 
derive a nonlinear system of the form (fT5b from a non- 
linear open quantum harmonic oscillator defined by L 
and 77, we proceed as in Section II to get the following 
definition. 

Definition 4.1 A nonlinear quantum system of the form 
( 1751 ) is said to be physically realizable if it is a repre- 
sentation of a nonlinear open quantum harmonic oscil- 
lator defined by a nonlinear coupling operator L and 
a nonlinear Hamiltonian H, i.e., if there exist H and 
L such that the matrices A(a,a i ),B(a,a^),C(a,a^) and 
Z)(fl,fl^) satisfy the following equations 

A{a,a < ) = i [V ,a]lL + i[H,a]\ 

/j(a,fl t ) = [l\a]l\ 

C(a,a^) = L; 

5(a,a f ) = 7; (23) 
where I — l. ^ ' 



The following theorem provides necessary and suffi- 
cient conditions for a nonlinear quantum system of the 
form ( fT5l ) to be physically realizable. 

Theorem 4.1 A nonlinear quantum system of the form 

( 1751 ) is physically realizable with T„ = I = ^ ^ 

if and only if 

[A(a,fl t ),fl t ] + [ajA^a 1 ") 1 "] 
+S(a,a + )7/j(a,a t ) t =0; 

[S(a,a t ),a t ] =0; 

[a,5(a,a t ) t ] =0; 

B{a,a^) = [C(a,a' ) T ,a] 7 and 

D(a,a r )=I. (24) 

In this case, the corresponding nonlinear self-adjoint 
Hamiltonian H is given by 

77 = (a t (r 1 A(a,a t ) - A(a,a t ) t 0" 1 a) (25) 

where h is as defined in Definition \2.1\ and the corre- 
sponding coupling operator is L — C(a,af). 

Note: Note that the physical realizability conditions as 
illustrated by conditions (l24l i impose the restriction of 
requiring C(a) to be a vector of first order polynomi- 
als of the annihilation operator a; i.e, the index m^, 

n 

in Definition 12. II is only 1 and C v (a) = £ $\p a p f° r 

v = 1 , . . . , m. Consequently, this in turn imposes the re- 
striction of requiring the matrix 7? (a , a ^ ) to be a constant 
matrix. 

5. Lossless Property of The Class of Non- 
linear Quantum Systems 

In this section, we consider a class of nonlinear 
quantum systems of the form (TT3T > and we give a quan- 
tum version of the lossless property for this class. 

Definition 5.1 Given an operator valued quadratic 
form 

r(P w j)=fQy + 2y f Sp w + ptRp w (26) 

where Q G C' 5x ' 5 , S G C^ x '" andR G W IX " 1 are constant 
matrices with Q and R Hermitian, we say that the sys- 
tem ( 1751 ) is dissipative with supply rate r([} w ,y) if there 
exists a storage function (a) which is polynomial op- 
erator valued mapping with 0(a) > Ofor all t such that 
0(0) =0 and 

/%(&(*) j(s))d S >$(a) (27) 
Jo 

for all t > and all operator valued j3 w (f ). 



Definition 5.2 The nonlinear quantum system of the 
form ( 1751 ) is said to be lossless if it is dissipative with 
Q = ~I, 5 = and R = I and 



r(p w (s),y(s))ds = 



(28) 



whenever 0(0) = <j)(a(t)) = for all operator valued 
p w (t) and all t > 0. 

Theorem 5.1 A necessary and sufficient condition for 
the nonlinear quantum system ( 1751 ) to be lossless with 
respect to the supply rate ( 1261 ) is that there exists a 
storage function 0(a) with 0(a) > for all t such that 
<j>(0) = 0and 



V0(a) T A(a,a T ) 



f\ - 



1 



B(fl,a t ) t V0(a) 



7-D(fl,a t ) + D(a,fl T ) = 



f\ - 



—C(a,a^) i( C(a,a' i ) 
— C(a,a^) 

0. (29) 



6. Physical Realizability And The Lossless 
Property 

In this section, we show that the nonlinear quantum 
system of the form ( TTBT l is physically realizable if and 
only if it is lossless with a specific storage function. 

Theorem 6.1 The nonlinear quantum system ( 1751 ) is 
physically realizable if and only if it is lossless with a 
differentiable storage function 0(a) satisfying 



[V0(< 



27 

-27 



(30) 



7. Example 

We consider a nonlinear quantum system described 



by 



A(a,a' ) 

B(a,a*) 

C(a,a ) 
D(a,a ) 



-k ai ai + 2a\a\ 



-ka^a^ — 2a^a ^ 




-J 2k, 



a 2 



/ 2k^a i 
f2kZa 2 



where a is the annihilation operator, a^ is the creation 
operator and k ai and k ai are constants. 



Hence, 

A(a,a T ) — 

B(a,a T ) = 

C(a,a T ) = 
D(a,a f ) 



—k ai ai + 2a\a\ 



ka\&\ 
~k(_i^ a 



- 2a\a\ 
-2a\ai? 
2 - 2a 2 a\* 



\/2kai 





\ /2k ai a i 
■sj 2k a J a j 



o 
o 






[2k n 







/2k, 



"2 



= 7. 



[B{a,a^),a T ] = 0, 
- [a,A(a,a t ) 7 '] = 



The matrices A (a,a y ), B (a,a y ) and C(a,a T ) satisfy 
the properties of the class of nonlinear quantum systems 
under consideration as given in Definition 12.11 In fact, 
C(a),a T ] =0and \A{a,a^),a T } = 
-2a\ 
2a\ 

Moreover, A (a, a^) and C(a,a^) satisfy conditions 
( TTol l in Definition ^. II In this case, h = 4. 

Also, condition (fTTT i in Definition 12.11 is satis- 



fied. In fact, 



2/7 



A{a,a^)^ 9 l a,a 



a l a 2 

* 9 

— ajflj 

9* 



and 



2^ [a^O 1 A(a,a J <) 1 a\ = 



— u ^ U2 

* 9 

aja\ 

9* 

— «2 "1 



hand, A(a,a^) — ^B{a,a^)C{a,a' ! ) 



Hence, 
1 

2fi 



On the other 



2fljfl2 

— 202^1 

2fl2*«l 

— 2fl|*a2 



A(a,a^0 'a, a \a^6 l A(a,a?),a\ 

i 2n 1 J 

A(a,a^) — — B{a,a^)C{a,a ! ). 



In addition, the system under consideration is physi- 
cally realizable. In fact, B{a 1 a^) = [C(a,a^,a] 7 = 



I2k„ 















-y/2K, 

-y/2k^ 

-^/2k^ 

On the other hand, [/i(a,a^),a^ 
a,7?(a,a t ) t ] =0 and 73(0,0^ =7. 



0; 



Also, 

-\a\a\ 

2* 
2 





2a\ 



4a*,a2 



-2a]* 



-2a\ 


Aa2a\ 

\a\a\ 

-2a 2 2 * 




[A(a,a^), 




and 



2a] 
-\a\a\ 

—4a*a2 



2af 



2a\ 


-4a2a* 





-2a\ 
\ci\a\ 



B(a 


a*)IB{a 


flt)t 





























- 



















Hence, [A(a,a^),a^] + [a,A(a,a''')'''] 

fi(a,a t )/fi(a,a t ) t =0. 

In this case, the corresponding nonlinear self- 
adjoint Hamiltonian H is given by 



H = -^(d Jf 9- 1 A(a,a r )-A{a,a J( )' i e- 1 ' 



2n 

■ 2* 2 

ia, a 9 - 



,2* 2 



and the corresponding coupling operator is L 
C(a,a T ). 

Moreover, the system is lossless. In fact, let 0(a) 

2a i 

2a 2 
2a\ 
2a% 



0(a) t = 2a\a\ +2a*,a 2 > 0. Hence, V0(a) 
Hence, [V0(a),a r ] 



2/ 

-21 

Also, V0(a) T A(a,a T ) = — C(a,a + ) + C(a,a' ) = 
—2k ai (a\a* + a[a\) — 2k ai (a2a\ + 0*^02) ■ 

Moreover, jB(a,a T )^V^(a) = — C(a,a') and / — 
5(a,a t ) + D(a,a t ) =0. 

Hence, the system is lossless using Theorem l5.ll 

8. Conclusion 

In this paper, the issue of physical realizability for 
a class of nonlinear quantum systems is developed and 
was found to be connected to the lossless properties of 
the system. This provides a connection between the 
quantum mechanical and nonlinear system theory prop- 
erties of the system. 
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